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ABSTRACT 

Let  U and  V be  independent  random  variables , and  let 
W = UV.  This  paper  concerns  the  distribution  of  U,  given 
that  V and  W are  distributed  according  to  the  gamma  distri- 
butions. It  is  shown  that  U is  distributed  according  to  a 
beta  distribution  if  the  distributions  of  V and  W are  central 
gamma  and  that  the  distribution  of  U is  degenerate  at  u = 1 
if  the  distributions  of  V and  W are  non-central  gamma.  The 
given  result  is  applied  to  determine  the  distribution  of  U 
when  V and  W are  normally  distributed. 
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1.  Main  results.  Let 


g_(x) 

am 


,m-l 

1ST  e ' 


x > 0 


denote  the  gamma  density  function  with  m degrees  of  freedom. 

The  non-central  gamma  distribiton  G . with  m degrees  of  free- 

m,  o 

dom  and  non-centrality  parameter  equal  to  6 ( > 0 ) is  given  by 
the  density  function 


(l.D 


V«(x)  * e 


*sr 

r=0 


*m+r 


(x)  It 


m-1 

- e"5'x(r)_2_ 


where  I (x)  denotes  the  modified  Bessel  function  with  para- 
m 

meter  m. 

Let  U and  V be  independent  random  variables  and  let 


(1.2)  W * UV. 

Suppose  that  the  distribution  of  V is  G . and  the  marginal  dis- 

Hlf  (J 

tribution  of  W is  G , , The  main  result  of  this  paper  concerns 

m f o 

the  distrubiton  of  U.  Clearly  U is  positive  with  probability  1. 
Moreover 


(1.3) 


P{0  < U < 1}  = 1. 


Otherwise,  let  P{U^l+£}  = a,  where  £ and  a are  positive 
numbers.  Let  c = (1  + £)~^.  Then 


Ee 


cw 


Ee 


cUV 


> aEe 


V 


(1.4) 


(2) 


where  E denotes  expectation.  The  left  hand  side  of  (1.4)  is 
finite,  whereas  the  right  hand  side  is  infinite.  Therefore 
(1.3)  is  true. 

It  is  easy  to  show  that  m'  <_  m and  6'  £ 6.  Let  E(U)  * 
E(W)/E(V)  = (m'  + 5')/(m  +6)  - y , say.  The  Laplace  trans- 
form of  the  gamma  distribution  is  given  by 

f0e'xx  dGm,5(x»  - (1  + *>~m  eitp<-rfr>-  x > -1- 

Therefore,  (1.2)  yields 


(1.5) 


(1  + X)“m’  exp  (-5'  + 15^x)  = E(1  + XU)"m  exp  (-6  + ± xu) 


>_  (1  + \y)  m exp  (-6  + j-  + x^.) 

by  Jensen's  inequality.  Comparing  the  two  sides  of  (1.5)  for 
large  values  of  X we  find  that  m'  < m. 


tla.bjx)  - 1 + §x  + bib  t 1)  It  + ••• 


denote  the  confluent  hypergeometric  function.  We  have 

EWr  " ‘’nm'f  r)  e‘6'  Mm'  + r,  m’  ; 5')  , r > -m’ 


Similarly 


r (m  + r)  -6  .v 

P (jn)  ® ( ^ + r , m,  o ) , 


_•  > -m. 


Therefore,  (1.2)  yields 


(3) 


(1.6)  EUr  = r^m  f(in*^  e~  5 ifni'+r  ,m' ; 6' ) / ^ r~~  e_<5,|>(m+r,m;6 ) 


= (6') 


3 m'  m3 

■'*  " T"  (6)  7 ~ T r (m*  -m)  /2  exp  (6j«l  _ 2/F 


(/6  - /T7)  ) (1  + 0 (r-1)  ) 

for  large  values  of  r.  The  asymptotic  expression  given  above 

is  derived  from  Formula  6.13.2  (12)  of  Erdelyi  (1953).  Since 

EUr  _<  1 for  r 2 0»  it  follows  from  (1.6)  that  6'  < 6. 

Let  H(u)  denote  the  distribution  function  of  U,  and  let 

f *(v)  = (1  + v)_m  and 
m 


f (v)  = 
m 


( v + u)  m dH  (u)  , v > 0 . 


Putting  v = t in  (1.5)  we  get 


(1.7) 


vm'“m  , 5'  . -5 

r- exp  (_  r+~G)=  e 


(1  + V) 


m' 


0(v+u)-m  exp  (^“)  dH(u) 


_ .-5  r»  ( - 6v)  r T(m)  * (r) 

- e ^r=Q  r!  fTm+rT~  fm  (v) 


= r(m)  e-<SH„  . (2/vTd)  f (v) 
m-l  m 


_ d 


where  D = denotes  the  derivative  operator  with  respect  to  v, 
and 


oo  , 2/nr 
(-x  / 4) 


H (x)  = l'  ______ 

r=0  r(m+r+l) r ! 


, 


(4) 


Note  that  J (x)  = (£)  H (x)  represents  a Bessel  function, 
m t m 


Therefore 


(1.8)  (x  -^-5-  + (2m  - l)-£-  + x)  H . (x)  = 0. 

cbc  ax  m_i 

Writing  the  left  hand  side  of  (1.7)  in  the  same  form  as  the 
right  hand  side  we  get 


(1.9) 


r(m,)vm’‘m  e"6'  H i . (2/eTrD)  f ,*(v)  = 

m “l  m 


T(m)  e”°  (2/06D)  fm(v) 


where  9 = v. 


First  let  m'  = m.  A transformation  of  (1.8)  gives 


(1.10)  ( 4x  + 4m  + c2)  H . (c/x)  = 0 

dx2  dx  m"1 

where  c is  a constant.  An  application  of  the  differential 
equation  (1.10)  to  both  sides  of  (1.9)  with  c = 2/e"6  gives 


49(6  - 6')  e~°  H . (2/91^0)  f (v)  = 0 

m-i  m 


(6  - 6')  fm  (v)  exp  (-p-^)  = 0. 

Therefore,  6=6'.  Hence,  P{U  = 1}  = 1. 

Next,  let  m > m'.  If  6 = 0 then  6'  =0,  since  6'  £ 6, 
as  shown  above.  Then  (1.6)  reduces  to 


_ T (m'  + r)  , r (m  + r) 

EU  'T(m')  ' r(m> 


fV 


dH* (u) 


(5) 


where  H*  denotes  the  beta  distribution  6 (u;m' ,m-m' ) . Hence, 
H(u)  = H*(u).  Suppose  that  6 > 0.  Writing 

e-6’  j. ^ (m-m-+r-l,  r(m+r)  Hm+r_^  £J+r(v) 

for  the  left  hand  side  of  (1.7)  and  applying  the  differential 
equation  (1.10)  to  both  sides  with  c = 2/05’  we  get 

— 6 1 r°°  ,m-m 1 +r-l.  _ , . ...  . d 
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r=0 


( r * ) r(m+r)  (40(5-6')-4r  ^) 


Hm+r-l(2/^)fm+r(v)  = 0 


or 


m'-m 


40(5-6') 


(1  + v) 


H-r  exp  <t4-7>  ♦ «««•  r r 


r=0 
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« (06  ' )s  (1+v) 


, m-r-s 


s=0 


s!  m+r+s 


= 0. 


The  above  equation  implies  that  6 = 6 ' = 0 , contrary  to  the  assump- 
tion that  6 > 0. 

The  foregoing  results  are  summarized  in  the  following  theorem. 

Theorem  1.  Let  W - UV  where  U and  V are  independent  random 
variables,  and  let  V and  W be  distributed  according  to  the  gamma 
distributions  ^ and  Gm,g1,  respectively.  Then  m'^m,  5'<  5 and 
P{0  < U <_  1}  = 1.  If  m'  = m then  5'  = 6 and  P{U  = 1}  = 1.  If 
m'  < m then  5'  = 6 = 0 and  U is  distributed  according  to  the 
beta  distribution  8 (u;m' ,m-m' ) . 


(6) 


Since  the  square  of  a normal  random  variable  is  dis- 
tributed according  to  the  gamma  distribution  with  a scale  factor 
we  obtain  the  following  corollary  from  the  above  theorem. 


Corollary  1.  Let  W = UV  where  U and  V are  independent  random 
variables,  and  let  V and  W be  normally  distributed  with  unit 
variance  and  means  equal  to  y and  y',  respectively.  Then 
y2  * y'2  and  P{U2  = 1}  = 1.  If  y = -(+)y'  / 0 then 
P{U  = -(+) 1}  - 1. 

An  extension  of  Corollary  1 is  given  as  follows:  Let 
a.,..., a , Z Z be  independent  random  variables  and  let 

x r x r 

Z = lP  a . Z . . Let  EZ.  = E. . 
i=l  11  11 

Corollary  2.  If  the  random  variables  Z,Z^,...,Z  are  normally 

2 2 

distributed  then  P{a^  = c^}  = 1 when  ^ 0 and  P{a^  = c^}  = 1 
when  E.  = 0 for  each  i = l,...,p  where  c. ,...,c  are  certain 

l l p 

constants . 

Proof . Since  Z is  normally  distributed  it  follows  from  the 
reproductive  property  of  the  normal  distribution  (see  e.g.,  Lukacs 
and  Laha  (1964)  Lemma  5.1.1)  that  a^Z^  is  normally  distributed 
for  each  i = l,...,p.  The  conclusion  of  the  corollary  follows 
from  Corollary  1. 

Let  a-,  , • . • , a be  m-component  random  vectors  and  let  A denote 
the  matrix  whose  ith  column  vector  is  , i = l,...,p.  Let 
Z.,...,Z  be  p independent  normal  random  variables.  Let 

r 


(7) 


EZ^  * Cj_»  Z = (Z^,...,Zp)  ' and  Y = AZ . The  random  vector  Y is 
distributed  according  to  a multivariate  normal  distribution  if 

and  only  if  X ' Y = (X'a.)Z.  is  normally  distributed  for 

i=l  1 1 

every  non-null  vector  X.  The  following  result  follows  from 
Corollary  2 . 


Corollary  3.  Let  a^,...,cXp,  Z be  independent.  If  the  distri- 
bution of  Y is  multivariate  then  P{X'ou  = c^}  = 1 when  ¥ 0 
and  P{(X'a^)  = c^}  when  = 0,  i = l,...,p  for  each  non-null 

vector  X,  where  c^,...fc  are  certain  constants  depending  on  X. 

Corollary  3 is  related  to  the  following  result  due  to 
Kingman  and  Graybill  (1970).  Let  Y^,...,Y  be  independent  and 
identically  distributed  random  variables  and  let  A = (a^j)  be 

a pxp  random  matrix  which  is  orthogonal  with  probability  1 and 
rP 

E(£  a..)  ^ 0 f°r  some  i*  Let  Y = (Y,,...,Y  )'  and  Z = AY. 

j=l  3 P 

Then  the  components  of  Z are  independently  and  identically  dis- 
tributed according  to  the  standard  normal  distribution  if  and 
only  if  the  components  of  Y have  the  same  distribution. 

Theorem  2 below  gives  a characterization  of  the  gamma  dis- 
tribution. Let 

F(t  h t)  _ ?°°  r (a+r)  r (b+r)  r(c)  xr 

f (a,b,c,x)  - -r(31  f(b)  r-(c+r)  pr 

denote  the  hypergeometric  function,  and  let 


(1.11) 


<f>  (X)  = F (a,b;c;-X) 


f tb_1(l-t)c“b_1(l  + tX) -adt 
J n 


a,b,c  > 0,  c > b. 


(8) 


It  is  seen  that  4>(0)  = 1 and  that  $(A)  is  a completely  monotone 
function,  that  is,  (-1)  <p v (X)  >.0,  A > 0.  Therefore,  <j>(A) 
represents  the  Laplace  transform  of  a probability  distribution 
on  [0,®).  A distribution  on  [0,®),  whose  Laplace  transform  is 
given  by  (1.11)  will  be  called  inverse-hypergeometric.  If 
b = c then  4>(A)  = (1  + A)  ais  the  Laplace  transform  of  the 
gamma  distribution  with  a degrees  of  freedom. 

Theorem  2 . Let  W = UV  where  U and  V are  independent  random 
variables,  and  let  U be  distributed  according  to  a beta  distri- 
bution 8(u;p,q),  say.  Then  W is  distributed  according  to  a 
gamma  distribution  with  m degrees  of  freedom  if  and  only  if 
m < p and  the  distribution  V is  inverse-hypergeometric , given 
by  the  Laplace  transform  <j>(A)  = F (p+2  ;m;p; -A)  . 

Proof.  Suppose  that  the  Laplace  transform  of  the  distribution 
of  V is  given  by  <j>(A)  = F(p+q;m;p;-A)  . The  rth  moment  of  V is 
given  by 


(1.12) 


EVr  = r (p+g+r)  r (m+r)  r(p) 

r(p+q)  r(m)  r(p+r) 


Therefore , 


(1.13) 


EvF  = EUrEVr 


_ F (m+r) 

r (m) 


The  right  hand  side  of  (1.13)  represents  the  rth  moment  of  the 
gamma  distribution  with  m degrees  of  freedom.  Therefore,  the 
distribution  of  W is  gamma. 


(9) 


(10) 
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